manuscript No. 

(will be inserted by the editor) 



John Baber 

Scaled Correlations of Critical Points 
of Random Sections on Riemann 
Surfaces 



Received: date / Accepted: date 



Abstract In this paper we prove that as N goes to infinity, the scaling limit 
of the correlation between critical points zi and Z2 of random holomorphic 
sections of the A^-th power of a positive line bundle over a compact Riemann 
surface tends to 2/(37r^) for small \/N\zi — Z2\- The scaling limit is directly 
calculated using a general form of the Kac-Rice formula and formulas and 
theorems of Pavel Bleher, Bernard Shiffman, and Steve Zelditch. 
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1 Introduction 

This paper studies the behavior of the critical points of gaussian random 
holomorphic sections of the iV-th power of a holomorphic line bundle L on 
a Riemann surface ill as — >■ cx), as is studied in [8], [9], and [10]. In the 
particular case where L — 0(1), the so-called hyperplane section bundle over 
M = CP^ , sections of correspond to homogeneous polynomials of degree 
N, the SU2 polynomials, so the results in this paper apply to the critical 

points of random polynomials ^ (^)cfc^'^ with Ck identically distributed 
gaussian random variables. In this way, this paper examines one small facet 
of the theory of random polynomials and random holomorphic functions. 

Since what may have been the first study of critical points of random 
curves in [15], this area of research has led to results of interest in mathe- 
matics, probability theory, and physics. For instance, the classical result of 
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Hammersley in [12] that for f{z) := X^j^g CjZ^ with Cj independent standard 
gaussian random variables, as A'' — >■ oo, the complex zeroes tend toward the 
unit circle in C and its generalization by Bloom and Shiffman in [6] (also 
discussed in [5]), namely that as -/V — >■ oo, the common zeroes of m random 
polynomials fk{z) :~ '^\j\<n '^^■^i^ ' ' ' '^"^ ^'^'^ concentrated near the 
"distinguished boundary" of the m-dimensional polydisc. Since the zeroes of 
a collection of m polynomials in m variables is almost surely discrete, for 
random fi, the set {,fi{z) ~ J2{z) — ■ ■ ■ — fm{z) = 0} is a random point 
process on C" of interest in probability theory. 

How much should zeroes and critical points of random polynomials or 
random holomorphic functions Xlfc°=o Id'^^ ^'^ expected to vary from their 
expected behavior? This type of question is addressed in [18], [19], [21], [20]. 
This paper examines how pairs of critical points are correlated by examining 
the 2-point correlation function, K2{z,w). 

The main theorem of this paper says that the scaling limit of the corre- 
lation between critical points of random holomorphic sections of the 7V-th 
power of a positive line bundle over a compact Riemann surface tends to 

as iV — >■ oo for small r := — C2|- i-S- 

Theorem 1 For any positive hermitian line bundle L over any compact Rie- 
mann surface M 

where r := dist(Ci,C2) i-e. as the distance between critical points gets smaller, 

2 

their scaled limit correlation approaches — 7; uniformly in Ci , C2 • 

K2ii{z, w) is calculated via the generalized form of the Kac-Rice formula 
of [14], [16] 

K{t) = J |C|JPD(0,C;t)rfC (2) 

whore JPD(a;,^; t) denotes the joint probability distribution of a; = f{t) and 
e = fit). 

Though we know no immediate interpretation of the constant 3;^, the 
fact that it is not is interesting. This contrasts with the fact that the scaling 
limit correlation of zeroes of random sections on a compact Riemann surface 
is O(r^) as was proved in general in [2] and [3] and specfically for CP^ in 
[13]. 

This paper is based on the thesis submitted to the Department of Math- 
ematics at Johns Hopkins University in 2010 which was read by Bernard 
Shiffman (Advisor) and Steve Zelditch. 

The introductions of [8], [9], and [3] give a description of the basic ob- 
jects of study and the physical motivation for them. The next few sections 
summarize the more thorough descriptions given there. 
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2 Notation and Formulas 

Throughout these definitions, M will denote a complex manifold of complex 
dimension n with complex coordinates {zi, . . . , Zn)- M will also be thought of 

as a 2n-dimensional real manifold with coordinates (xi, yi, . . . , .t„, y„) where 
Zj =: Xj + iyj. C denotes Lcbcsguc measure on C and B denotes the borel 
subsets of C. In general L will be a holomorphic line bundle over M. For 
standard results and definitions about line bundles, see Chapter 1 of [11] for 
instance. 

Wo begin by summarizing our notation. Tm denotes the set of smooth 
complex-valued vectors on M. i.e. Tm,p is the space of C-linear derivations 
in the ring of complex-valued C°° functions on M near p. 



Tm = T'm ® T'. 



' M 



(3) 



where 
T'm ■= span 



d 


d 


dz\ ' " 


dZn 


d 


d 


dz\ ' 


' dZn 



the "holomorphic" tangent space (4) 



Tj^ denotes the dual space of Tm, i.e. the set of smooth complex- valued 
covectors or 1-forms on M . 



Tm' 



M 



where 

Tm '■= span {dzi, . 
Tm ■■= span((izi, . 

Recall that 
where 



, dZn) 
, dZn) 



the "holomorphic" cotangent space (6) 
the "antiholomorphic" cotangent space (7) 



d = d + d 



i=i 



(8) 
(9) 



df 



df 



We let AP{L) denote the sheaf of smooth L- valued p- forms, i.e. for any 
open UcM, Ap{L){U) := 



{u;\^^eu{x) I we /\tX^ 



and eu a local frame above U}. (10) 



We also let AP''^{L) denotes the sheaf of smooth L- valued (p, Q')-forms. i.e. 
for any open U cM, AP'1{L){U) := 



{uj\®eu{x) I we T^P'^^ 



and eu a local frame above U} (11) 
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The d operator is extended to act on sections via 

w e I — > 9w e 

Using the T]^ = T^' decomposition, we write V = V + V" for any 
connection V where 

V (12) 
W" : A\L) ^ A'^'^L) (13) 

Given a hermitian line bundle (L, ( , )h) — >■ M, the Chcrn connection 

associated to ( , )h, will be written V/j. When the h is obvious, V?i and 

{ , )h will just be written V and ( , ). 

For given N , we will choose local coordinates on U and a local frame 
for over U C M such that h{z) = 1 - [zp + 0(|z|3) by taking an 

arbitrary frame and multiplying by a smooth function with appropriate first 

and second order terms. Then 



and 



=d + Ndlogh (14) 
|^logft=|jlog(l-|.|^ + 0(|.H) 

= il + 0{\z\'))i-z + 0{\z\')) 
= -z + 0{\z\') 

Also 

K = d (16) 

The curvature form for will be written Oh- Note that because i is a 
line bundle, Oh is just the 1x1 matrix \dO\ = \dd log /i] . 

We now summarize notation from probability theory. For a random vari- 
able 

X ■ (I?, P) (E, BorclSets, dLcbcsgue) (17) 
we'll write the cumulative distribution function of X as 

Fx{t):=V[X-\{-^,t])] (18) 
When X e L^{P) we denote the expected value of X by 

E[X] := f XdP (19) 

Note when X happens to have a probability density function fx, 

E[X]= [ tfx{t)dt (20) 
Jm. 
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Definition 1 A centered complex gaussian random variable is a random vari- 
able 

X : (r2,r,P) (C,S,£) (21) 

whose distribution is 

iX,P){B):=P[X-\B)]= [ -l-e-^\'\" dC{z) (22) 

Jb ■^o- 

When CT = 1 we say X is a standard complex gaussian. 
Note any centered gaussian has expected value 0: 

/ XdP^ [ ^e-^^'^' dC{z) = (23) 
Jn Jc '^'^ 

Definition 2 More generally a collection of random variables Xj : f2 ^ C 
is said to be jointly gaussian if the complex valued random variable 

aiXi + ■■■+ anXn (24) 

is a centered complex gaussian for any aj £ C. 

Definition 3 The n x n symmetric positive semi-definite matrix 

A:=[E[Xaj]],,^,,„„ (25) 

is called the covariance matrix of X. When the Xi are linearly independent, 
as in our calculation, A is non-singular, i.e. positive definite. 

When the Xi are linearly independent, Definition 2 is equivalent to a 
more probability density style description. Specifically, Definition 2 in this 
case is equivalent to demanding that the random vector 



: f2" ^ C" (26) 



has distribution 

(X.P)(B) = P[X-\B)] = _2_e-<^-^-^> dz (27) 

where 

dz = d/:{zi) A . . . A dC{zn) (28) 

The distribution X«P for any P as above is called the joint probability 

distribution of the Xj . 

Lemma 1 If Xi,. . . ,X„ are jointly gaussian then the entries of L{X.) are, 
too, for any linear surjection i : C" -» C"* 
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Proof If [L] = [iij] and 







\ 




V 




J 


Y 



then aiYi H h amYm 

V + • • • + ^i„X„) + • • • + am{^mlXi +■■■+ i^nXn) 
= {ai£ii H h am^ml)-'^! -I 1- {aihn H h am^mn)-'^n 

a centered complex gaussian for any aj e C. □ 



3 Random Sections cind the Two Point Kernel 

Here we define what we mean by "random sections" of the bundle M. 

Definition 4 The metric h induces hermitian metrics /i^ on given by 

h'^{z) := h{z)^ i.e. (si • • • (g) sjv, ti <?> • • • iAr)z,N = 

{heu (g) • • • fNeu,gieu • • • gNeu)L^ 

=/l • • • /nSi ■ ■ ■ gN{eu • • • e;/, 6(7 • • • eu)L'^ 

:=/i • • • InSi ■ ■ ■ 9N{eu, eu)L = /i • • • InQi ■ ■ ■ 9Nh{z)^ 

Definition 5 Using we can create a new hermitian inner-product on 

H°{M,L^) by 

{s,t)^[ {s,t)LN(No\M s,tG H°{M,L'^) (30) 
Jm 

Throughout the rest of this section, (s^) will denote an orthonormal basis 
for H°{M,L^) with respect to ( , ). 

Definition 6 We can define a gaussian probability measure P on 
H°{M,L^). Given 

e 

H%M,L^)5s = Y,Cj{s)sf (31) 

for any borel collection of sections S, 

P[S]:= J l-e-^"'^')'"^^')) dc{s) (32) 
where dc{s) is 2i!-dimensional Lebesgue measure. 
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P is characterized by the property that the 2i real variables Re(cj) and 
Im(cj) arc independent random variables with mean and variance 1/2. 
Specifically 

E [cj] =0 E [cjCk] = E [cjCk] = Sjk (33) 
For ci , . . . , Q jointly gaussian, consider the random holomorphic section 



s{z) :=Y,Cjsfiz) 



(34) 



and the map 



Cl 



■ (v;.)(z) ■ 
(v;v;s)(z) 
(v/v;s)(z) 

(V'^s)(ti;) 

(v:„v:„s)(ti;) 

.(VjV:„s)(«;). 



E;=ic,(v;.f)(z) 
E-=ic,(v;v;.f)(z) 
E}=ic,(v;'v;sf)(z) 
E5=ic,-(v;sf)(t«) 

E-=ic,(v;v:„.f)(^«) 
E5=ic,(vjv:„.f)(^«)j 



(35) 



For fixed z and tz;, A is a linear map so Lemma 1 says the entries of A(c) are 
jointly gaussian. Their covariance matrix is 



where 



and 



^=[A-^]L=i (36) 
= E [Xj O (37) 

Ai = (V»(z) 
A2 = iVXs)iz) 

A3 = i^:V:s){z) (38) 
A4 = {Ks){w) 
A5 = (V;V'^s)(w;) 
Ae = (VjV;s)(«;) 

by abuse of notation. Each entry should be replaced by its coefficient when 
written in a local frame about z and w. 

In fact all of the entries of A used in our calculation can be rewritten 
in terms of derivatives of an important invariant of P called the "two point 
kernel" . 

Definition 7 The two-point kernel (or covariance kemet) for H°{M,L^) is 
defined by 



nN{z,w) := ^sf (z) (g) s^{w) G if > 



{z,w)€M 



(39) 
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Since is hermitian, 11 n is the Szcgo kernel of (i^, h^), i.e. the orthogonal 
projection 

nN,h-,VoXM ■■ L\M, L^) ^ H\M, L^) (40) 
with respect to ( , ). 

Tliv and the entries of A are related because 



E 



s{z) ® s{w) 



= E 



^c,sf(z)®^Cfesf(«;) 



i=i 



^ E[c,cfc]sf(^)®sf(«;) 



j,k=i 



= 5,,sf{z) ® s^{w) = ^ (^) ® sf («;) 

=: nN{z,w) 
so differentiating both sides yields, for instance, 



y^V^HNiz, w) = V;VJ,'E s{z) s(u;) 

= v;e 

= Ai4 



= v,'e 



.s(z) (g) VJs(i(;) 



5(2:) (g) V^„s(u;) =E \7's(2;) ig) V'^s(w) 



(41) 



(42) 



4 The Kac-Rice Theorem 

Various generalizations of Rice's original theorem [15] (3) are referred to as 
"The Kac-Rice Theorem" in current literature. What is meant by "using the 
Kac-Rice theorem" is that the expected density of zeroes of a random linear 
combination of functions 



/a (a;) :=^aj/j(a;) 



(43) 



is found by integrating the joint distribution of /a and with /a replaced 
by against ||/al|. For instance, the single- variable real Kac-Rice theorem 
says the following: 

Take P a probability measure on and fi,---,fe a collection of analytic 
functions. For fixed t, 

Xt := /a(i) : ^ K (44) 

and 

6 := m : ^ K (45) 

are random variables so they have a joint probability distribution, namely 
the distribution for the random variable 



Xt := 



Xt 

it 



(46) 
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specifically 

{Xt),P=:Dt{x,0 (47) 
Definition 8 For the function define the measure 

/afe)=0 

The demand that the fj be analytic ensures /a has only finitely many zeroes 
on bounded intervals, so Zf^ is a sum of point masses. Zf^ can be gener- 
alized to mean the current of integration along the regular points of the 
variety {/a = 0} even when not discrete, but that is not necessary for our 
computation. 

The Kac-Rice theorem says 

E[ZfJ = K{t)dt (49) 

where 

K{t):= f D,{0,mdC (50) 
Jr 

Definition 9 The K(t) dt above is called the one-point correlation or one- 
point density oi Zf^. 

Definition 10 If we define the measure of n simultaneous zeroes 



E 

{teMx-xM|/,(ti) = -=/a(t„)=0} 



Zl--= E '^t (51) 



the Kac-Rice theorem says the same for a measure Kn{t) dt called the n-point 
correlation or n-point density of Z^^ . 

This paper is, in fact, concerned with the 2-point correlation of the si- 
multaneous zeroes of two random sections Vs^{z), Vs^{w) £ (g) for 
s^ e H^{M, L^) which are called critical points of . Here "random" means 
that the arc chosen with respect to the gaussian probability measure on 
H^{M, L^) given in Definition 6 and dime M = 1. In this particular case, 
the Kac-Rice theorem is 



E [Zvs(z),Vs(w)] = K2{z, w) dz dw (52) 



where 



K2iz,w)= D,,^{o,odet{eier)^det{eier)^d^ (53) 

Jw 

= [ £'.,^(0,0 • I det^^l • I det(C')| d^ (54) 
Jw 
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due to the makeup of and Here ^^^^(x, ^) is the joint probabihty 
distribution of the random sections 



and 



e 



xi = Vs{z) = Vls{z) 
xi = Vs{w) = V^s(w) 

a = vvs{z) 

6 = V"Vs(z) 
6 = VVs{w) 
^4 = V"Vs(w) 



since s is holomorphic 
since s is holomorphic 



^'VVs(z)" = 



V'Vs(0) V"Vs(2) 



V"Vs(^) V'Vs(z) 



V'Vs{w) V"V.s(w) 
V"Vs(w;) V'VsH 



V/V>(^) V/V>(z) 



v,"v>(^) v;v;s(z) 



(55) 



as in [4] (33). Here 
W = 



-N\ 



M 



-N\ 



xV{{n,^L^)JxV'{{T*^^L^)J 



and means Lebesgue measure with respect to the hermitian metric on W. 

As in the previous section, xi, 0:2, ^1,^2, ^3, and ^4 are jointly gaussian 
with covariance matrix A =(37) so 



1 



TT^ detZ\(2:,u') 



exp 



so 



£'.,«,(0, 0,^1,6,6,^4) 

Dividing A into blocks 
A 



TT^ det A{z, w) 



exp 



A-\z,w) 
A-\z,w) 



x 




X 










) 













c 






) 



^2X2 [^]2X4 
[-^14x2 [<^]4x4. 



(56) 



6x6 



and using the formula for inverting matrices presented in blocks 

I -A-'^B] [A-i 
/ 





I 


0" 




-B*A- 


1 / 



A-^ + A-^BA-^B*A-^ 
-A-^B*A-'^ 



2x2 



- A-^BA-^ 



4x2 



2x4 



4x4 



(57) 
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meaning 






2x1 



4x1 



A-^ + A-^BA-^B*A-^ 
- A-^B*A-^ 



-A-^BA-^( 



A-^BA-^ 
A-^ 



so 



1 
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) 



[{-A-^BA-^^y {A-^0*] 



Using the formula for determinants of matrices presented as blocks says 

det^ = (det^)(detyl) (58) 
so the integrand in K2{z, w) is 

^-{A-\z.w)i.O 



TT^ det A{z, w) det A{z, w) 



(59) 



Now, as mentioned above, the particular entries of the ^ allow the 
det(^^(^^)*)5 to be simplified as in [8](34). 

= det 



( 


'a 6' 




"a 6' 


^ dot ^ 


"a 6' 




"a a' 


■) 




a a. 




6 a. 


6 a. 




a a. 





(lap-iany ((laP-iaP)')^ 
lap-ian-liap-ian 
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Identifying W with C'', the two-point correlation of critical points on a Rie- 
mann surface is 



5 The Scaling Limit 

Although the critical point equation Vs{z) = is not holomorphic, it is still 
smooth, so the results of [3] about the zeroes of random smooth sections 
apply. The main theorem (3.6) of [3] actually says that as iV ^ oo, the 
"scaling limit" of the correlation of zeroes is independent of choice of M, L, 
and h. Specifically 



(61) 



where K^j,^^^ depends only on n, k, and m. 

If we write the 2-point correlation (60) as K2 to reflect the N dependency 
in our case, this theorem says 

so proving Theorem 1 only requires that we calculate^ the; N limit of the left 
hand side for z close to w for a particularly nice choice of M, L, and h. 

Theorem 3.1 of [3] says roughly that in the N limit, the 
■^]Jn{^^, ■^=) entries in can be replaced by II^{z, w), the Szcgo kernel 
of the reduced Heisenberg Group Hj-ed which we define below. For a more thor- 
ough geometric discussion of Hred's construction and properties, see [3]§1.3.2. 

Definition 11 Take the trivial bundle L := C x C over C with curvature 
h{z) := e~l^l . Then h~^{z) = e'^' gives a metric on the dual bundle L* — ^ C. 
Form the "circle bundle", X, of elements v G L* such that h~^{v) = 1. i.e. 

X = {{z,C) e C X C I Id = e-T-}. This bundle X ^ C is the reduced 
Heisenberg Group, written Hrod- When necessary, since X = C x S'-'^, we will 
write elements as {z,6). Because L — )■ C here is the trivial bundle, we may 
use the frame eu = ec = I, the constant function 1. 

The Szego kernel for Hj-cd: -^i^ is by definition the kernel of orthogonal 
projection from £^(Hred) to the Hardy space for Hredi ^i- These spaces can 
be viewed as 

£2(H,ed) = |/>,^) = f{z)e''e^ f e , jj{z)e-\^\' dzdz < ooj 

(63) 

and 

nl = lf{z,0)=f{z)e''e-- 



/ holomorphic, y f{z)e (iz(iz<ooj' 

(64) 
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In fact, (pk{z) := "7=^1^'^^*^^ ^ orthonormal basis for £^(Hred) since 

-J^lj,'e'<'-^>e-»'[t,.4.)=S,, (65) 

SO the kernel for the orthogonal projection £^(Hred) — >■ "Hf is 

°o 1 

n^iiz, 9), {w, ^)) = V <fk{z)n{w) = -e^(^-^)e^'^-^l^l'-^l-l' (66) 

k=0 

The extra factor e^l^l "ll""! appears because we have chosen a non-trivial 
metric h for the trivial bundle. As long as the connection is computed cor- 
rectly, we could work in any frame, so the formula becomes 

n^{z,w) = -6'"^ (67) 

TT 

based on the frame ec = 1 (see Definition 7). Since (60) doesn't change when 
nf^{z, w) is multiplied by a non-zero scalar (hence likewise if A is multiplied 
by a non-zero scalar), we will use 

n^{z,w) = e'^ (68) 

for ease of calculation. 

The Chern connection is defined by its action on a frame 

Vf ec(2:) = (V^)' eciz) = -zdz® ec{z) (69) 

Often, we will write n^{{z, 6), {w, ip)) and mean only the function coefficient 

For our case, we form the Szego kernel for the A^th power of an arbitrary 
positive line bundle over a Riemann surface {L^ , h^) — )• M similarly, defining 
the circle bundle 

Xm~ {seL* \{s{z),s{z)) = l] 
and calling the Szego kernel 

Hn-XxX^C (70) 
This Szego kernel is related to our earlier two-point function by 

J7„(-^,-^') , h{z)t{e},)"{z)®h{w)t{eir{«,) 



With those definitions, we are in the position to state Theorem 3.1 of [3] 
precisely in our case: 

Choose zq € M, local coordinate map z, and a local frame bl over a 
neighborhood of zq so that 

0h{zo) = iddlogh){zo) = dzAdzU, (71) 
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and 

dh , , d'^h , , „ 

^(.o) = ^(.o) = (72) 

then 

1 -- / / U d\ ( V if 

lif{{u,d\{vM)^o[-^ (73) 

where o[^^^ means a function whose C'^-norm is o{^-^^ in the standard 

sense for all k. 
For 

U{z) := IIn ( (zo + ^, 0) , (zo + w^, 0) ) (74) 



5«;(^) := (^^J = ^fjv(^(2:o + ^,0) , (^0 +t«,0) ) (75) 



[3]'s theorem says 

Taking derivatives on both sides of (75) 



!^[^\= (-) = ^r((-' 0), («;, 0)) + O ) (76) 



IN VviV 

so 

= |j(ilH((.,0),KO)) + o(-^)) (80) 
Since we chose h so that ^ log/i = —z + 0(|2:p) , 



iV ^ + o — 



|2 



-\/]V^ + 0(| 



z?) 
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= 0), 0)) + o (-^) ) - 1. + ^od.n 

= VfiT«((^,0),KO))+o(-^) 

-^K' ■■= ]^V;VjiT^, 

= (V?)' (V«)" 77r((., 0), {w, 0))|(,^,,^,) + O (;^) 

^^""-'^ ^V;VJVJ77^ (-^, |(C,.c,0 

= (V«)'(V«)"(V«)"iTP((.,0),KO))|(,^,,^,)+o(-^) 

= (V»)' 77r((., 0), O))|(c,,c,o + O 

T^^p'l ^= ]^v;v;vjvjiTA, (^-^, 

= (V«)' (VH)' (VH)" (V«)"77H((.,0),K0))|(,^,,^,) +o(-^) 
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w 



(V«)' (V?)' n^iiz, 0), 0))|(c^,c^,) + ^ 



1 ^P2 



A/-3 



7V3 



w 



N VN 



kcp.Cp') 



= (V^)" (V»)"i7H((.,0),K0))|(^,,f^,)+o(-^ 



1 /Of 2 ._ 1 rr 

]^V2 — ^ 



TV VN 



\/]V 



and 



where is A with all of the 77/^ terms replaced by 77^ terms and 



C = 



(81) 



So, finally, 
1 



7r6iVi6 detZ\H 



JV-=*((/l«)-i£,£> 



II detail (82) 



where is zi with all of the IIn terms replaced by iiP terms and || det^|| is 

shorthand for 

|det^i||dct^2|^ 

Now perform the change of variables 



Xi 
X2. 



:= 



i4. 



(83) 
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Now II detail = N^W dctvll and = N^'^dv so 

J-Ki^ (^,^)= ^ [ e-<(^")-^-'V>|| detvll rfv (84) 



6 Additional Definitions and Notes 

Definition 12 A critical point of s ^ H'^{M,L) with respect to V/i is any 
z e M such that Whs{z) = 0. 

For almost any s e H'^{M, L), the set 

Crit"^" (s) := {z e M I V/is(z) = 0} (85) 

is discrete, so the following definition makes sense. 

Definition 13 The measure associated to Crit^'' (s) is 

zeCrit'^'i(s) 

where 6z is the point-mass at z. 

Definition 14 The volume form associated to h, dVh, is given by 

dVh-^i-'-ddlogh) (87) 



m! V 2 

Definition 15 The two-point correlation was not originally defined as the 
necessary integrand in the Kac-Rice formula. It can be defined directly as 

, E [#[Crit^''(^) n B,{z)] . #[Crit^''(^) n B,{w)]] 
K2{z,w) = hm — ^, , rp . ^ D , (88) 

where 

B^{z) :=the ball of radius e about z 
#A :=the cardinality of A (89) 

K2{z,w) also comes from the distribution equation 

E [Cj'^ mcj'^] = K2{z,w) dVh{z)mdVh{w) (90) 

where M is the product on currents defined in [17] by 

SMT = ttIS Att;t gV'p+\M X M) (91) 

for S € V'P{M) and T € X>"?(M) where tti, ttz : M x M M are the 
projections to the first and second factors, respectively. 
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Definition 16 Throughout the calculation for n € [l,oo), will mean 

a function / such that 

3d,M>0 such that (^t G (0, 6) =^ \f{t)\ < Mt"^ 

In fact, every time it is used in this paper, it is sufficient to think of 0(t") as 
a function real analytic at whose first non-zero Taylor term when expanded 
there is a multiple of t". i.e. 

fit) = a„r + a„+ir+i + a„+2r+2 + • • • 

with a„ ^ 0. Technically, any time 0(t") is mentioned, it would be necessary 
to mention the radius of convergence, and often manipulation of a term 
involving 0(t") will result in a new term involving 0(f"), where the radius 
of convergence has shrunk. This calculation only requires that the radius of 
convergence stays positive. As long as this is the case, the reader should not 
pay attention to this technicality. 



7 Proof of the Main Result 

To prove Theorem 1, we want to find 



(92) 



for any (1,(2 & M where dime M = 1. As above, Kac-Rice says we need only 
calculate 

1 ('^]= —L— [ e-<(^")-^-'->|| detvll dv (93) 

' ttS det ZiH II II \ J 

As in (58), detzi" = (det A")(det yl") and since K^j^{z,w) depends only 
on the distance between z and w, we can choose z = and w = r > 0. So 



Kr,,{z,w) = Jir) := 
1 

776 det(^(0,r))det(yl(0,r)) Jc*' 
where 



/ 



\hi\'-\xi\ 



\h2\^-\x2\^\e-^^"^''''->'^Uv 



(94) 



The absolute value bars simplification using Wick's formula as in 
[3], the fact that Vs is not a holomorphic section bars using the Poincare- 

Lelong formula as in [2], and, unfortunately, we are unable to use the inge- 
nious method used in the proof of Lemma 3.1 of [9] where the authors were 
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able to rewrite J using Fourier transforms. In [9], the authors noticed that 
they could replace each — |a;j|^| by 



lim 



3 



(95) 



because it can be simplified to 



e'^^i^i dp 



lim — [ \p\e-''M" ( [ lim e*«(P-l'».l'+ 

= lim ^ / \p\e-'>^' ( [ lim^e*«(P-l''-'l'+l^^l''dn dp 
= hm / b|e-^ilfl'(27r5o(p - |/i,f + )) dp 

= lim / Iple'^'^IPl'^ift |2_| |2(p)cip 

= lim ||/i,f -k.f |e-^ill''^-|'-l^^-|'r 

e'-i-O ' ' 



\'-\xA'\ 



With that substitution and some work, 



J = hm Je = — J 



biP2|e 



«(ClPl)p«(?2P2) 



det(«ylL' - /) 



■ d^i d^2 dpi dp2 (96) 



with 



D 



-6 



(97) 



In the case of [9] , the authors took advantage of the fact that det{iAD — I) 
was the product of many linear factors and the integral could be done using 
residues. In our case det{iAD — I) is an extremely complicated rational func- 
tion of r and e'^ . So in this paper we will carefully expand J as a function 
of r. 

Calculate A as in [3] and [8]: 
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E 













A:=C - B*A-^B 
A := [Al,] 

B := [Bl,,\ = E 
C := [C- ] = E 



Al A\ 
Al Al 



^11 -^12 -^21 -^22 



d2 d2 d2 d2 
-''ll -''12 -''21 ^22 

-/ol 1 /ol 1 /ol 1 /ol 1- 

*-^21 '-^22 



/-<1 2 /^l 2 (^12 (^12 
"-^11 '-'12 '-'21 '-'22 



/^21 /-i2 1 ^2 1 ^2 1 
^11 '^12 '^21 '^22 



ri22 ri2 2 /^2 2 ri2 2 
y\l ^\2 ^21 ^22 . 



where p,p' G {1, 2} q, q' G {1, 2} and the p, q index the rows, and the p', q' 
index the columns. 



CP 1 
'^p'2 

CIA 

^p'2 



= v;\7;vj7jf(^,^i;)l(^^,g,) 

= V;<(z,«;)|(c,,c,0 

= V;v;\7:'Vji7f(2,«;)|(^^,g,) 

= V;v;7Jf(z,u;)|(c,,c,,) 
= y:X'7Tf(z,«;)|(^^,^^,) 



Notice p, q index rows and p' , q' index columns. 

For any function II{z,w) holomorphic in z and antiholomorphic in w 
where 

V^eiz) = V^e{z) = g{z) dz (g) e{z) 



VzHz) = V'_'e{z) = g{z) dz (g) e{z) 

V;'e(z)=0 

V;e(z) = 
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we have 

%{n{z, w) ® e{z) O e{w)) 

e{z) (g) e(w) + n{z, w) i 



an , 

oz 



g{z) dz (g) e{z) 



(g) e(w) 



— h g{z)n{z, w) ) (g) (g) e(2;) (g e(u') 



\7'\7'(77(z, {g e(z) e(u;)) 



ail 
dz 



+ g{z)n{z,w) 



g{z) dz (g e{z) 



e{z) e(w) 



(g e{w) 



IS) dz (E) dz (E) e{z) (E) e{w) 
__ + fn{z,w) + 2g(z)^ + H 
(g dz (g (i^; (g e(2:) (g e(t«) 



Vj(7I(2, u;) (g e(z) ® e{w)) 

e(z) (g e(w) + TI{z^ w) (g e(z) 1 



9n 

-—r aw 
ow 



g{w) dw (g e(w) 



dn 

+ g{w)n{z, w) ] E)dwE) e{z) (g e{w) 



VjVj(7T(z,«;)0e(z)0e(w;)) 
/57T 



(g dui 



e(z) (g e(w) 



+ g{'w)n{z, w) ) (g dtz; (g e(2:) 1 

dn 



g{'w) dw ® e{w) 



\ dw 



E> dw (E dw E) e{z) E) e{w) 
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XV;;,{n{z, w) ® e{z) ® e{w)) 



dn 



i^dzdw ^^^^ dz 
— + g{w)n{z,w) 



dn 



e{z) (g) e(w) 



g{z) dz (8) e{z) 



(8) e{w) 



\ dzdw dw dz 

® dz® diB® e{z) ® e{w) 



■ g{z)g{w)n{z,w) 



/ d^n dg dn -r-^ d^n 
ydzdw'^ dw dz dzdw 
®dw®dw® e{z) ® e{w) 
d^n dg ^. , ^-r^dn 



+ g{w) \ ®dz 



+ 



g_^n{z,w) + 2g{w)—+g{w) n{z,w] 



dw"^ 

® dw®dw ® 



g{z) dz e{z) 



e{w) 



d^n , ,d^n 

+ 9{z) 



dzdw"^ 



d^n 

2g{w)^^^ + 



dw"^ ' 'dzdw 
f dg_ 



dl 
dw 



+ g{w) 



+2g{z)g{w)^ + g{z) (^^+g{wf^ n{z,w) 



® dz ® dw ® dw ® e{z) ® e{w) 



dn \^ 

dz 

J 



In this particular case where g{z) = —z because h{z) := e n{z, w) 
e^^, (i := 0, and (2 '■= r, we have 



Ap, = e^™ {1 + zw — zz — WW + zw) 
B^i-y = e^""' (z — w) (zw — zz + 2 + zw — wiu) 



e^"' (w - z) 



2 — 4 zz — A WW + 4 zw — 2 wzz — 2w zw — 2 wzw 
-2z zw + Azw + w z +w w +z z +z w + Awzzw 



CP^,l=e^^{w-zf 
Cf^ = e^- {z - wf 

/op2 ziv 



(98) 
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so 



1 



so 



where 



B = 


r 

r (2 - 7-2) 


-r 
— r 


(2- 



r^) r 









2 





— 4r 


2+^4 


+ 2 




C = 




-4 + 


1 

+ 2r2 




^2 




1 






r2 


1 












Air) = 



1 - 2 r2 + 



iM2i(r)=r2 (-2 + r2) 

iMi2(r)=r2 (-2 + r2) 
iM22(r) = -e'^' + l-r2+r4 
iM32(r) = -r^ (e'-' + 1 - r^) 
iM42(r) = -e'"' + 1 - 

iMi3(r) = 4 r^e'"' - 4 + 3 - - r^e'^' 

'M23(r) 

'M33(r) 
'M43(r) 



2e'' +2 



-2e'' +2 



= r 



_^2 + 1 _ ,.2^ 

-e'-' (2e''' -2 + 2r^-r6) 
2 (_2 + r^) e'-' 

2 I , ^ „2 ' 



iM24(r) = -e"-' + 1 - 



(99) 



which is actually a function oit := r"^ > 0. 
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I.e. 



where 



2Mii(i) = -2 e* + 2 - 2 + 
^M2i{t)=t {-2 + t) 

2M3i(i) = 4i e* - 4i + - - i^e* - 2 e* + 2 
'^Mii{t)=t (-e*-l + t) 

2Mi2(t) =t (-2 + t) 
2M22(t) = -e* + l-i + i^ 
=t (-e*-l + t) 
2M42(i) = -e* + 1 - i 

2Mi3(t) = 4i e* - 4i + - - i^e* - 2e* + 2 
^M23{t) = t (-e*-l + t) 
2M33(i) = e* (-2 e* + 2 - 2t2 + 
'M43(t) = i (-2 + ^) e* 

2Mi4(0 = * (-e* - 1 + 

2M24(t) = -e* + l-t 
2M34(i) = t (-2 + t) e* 
2M44(t) = e* (-e* + 1-^ + ^2) 



so 



det A{t) = 

{e*ft^-eH^-4t^e*-4t^ {e*f+12t^ {e*f-12eV-12e* + 12{e*f-4{e*f+4 

-e* + 1 - 2t + ^2 

^ -t^ + J—t^+ t^" + ^^t"+ ^^^^ ti2^Q/^m (100) 

1 'isisja /inao/inn QORcioo ooQOQOsn ' V / v y 



6480 3888 4082400 326592 29393280 

det A{t) = e* - + 2t - 1 

= 3t-lt' + lt^ + ^/ + 0{t') (101) 

Now 
yl-i(t) = 

-t^e* +t4e2t - 4t3e* - 4t3e2* - 12f2e* + 12i2e2t - 4e^t_ i2e*+ Ue'^t+A ' ^^^^ 
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where 





[t) 




it) 




it) 




it) 




it) 


^M22 


it) 


'M32 


it) 


3M42 


it) 


'Mi3 


it) 


'M23 


it) 


^M33 


it) 


^M43 


it) 


^Mi4 


it) 


^M24 


it) 


3M34 


it) 


3M44 


it) 



SO 



where 



<''(e')'-t2(e*)^-e*t2+4t (e*)'-4te'-2(e*)V4(e*)^-2e* 

t'e' + f {e*Y + eH^ + 4te' -4t (e*)' + 2 (e*)'-4e* + 2 
-2*^6* + 2i (e*)^-4ie* + 2i 
-tV-t3 (e*)' + 2i2 (e*)'-2e*i2 

(e*)^*-e't*-4t3(e*)^-2 6**2 + 10 *"(e')'-4t(e')%4te'-4(e')^+8(e*)^-4e* 

-2t^e^ + 2t (e*)^-4ie* + 2t 

4 i^e* + 2 _ 10 + 4 1 e* - 4 i + 4 (e*) ^ - 8 e* + 4 

i3e* + <2 (e*)' + e*i2 + 4ie*-4t (e*)' + 2 (e*)'-4e* + 2 

-2*^6* + 2i (e*)^-4fe' + 2f 

t^eJ - e*<2 -t2 + 4ie*-4t-2 (e*)^ + 4 e* - 2 

-t^ - t^e* + 2 e*t^ - 2f 

-2t^e^ + 2t (e')^-4te* + 2i 

4 t^e* + 2 - 10 6**2 + 4 1 e* - 4 1 + 4 (e*) ^ - 8 e* + 4 

-f''-f''e* + 2 6**2-2*2 

e**'' - - 4 t^e* - 2 12 + 10 eV - 4 i e' + 4 1 - 4 (e*) V 8 e' - 4 



yl-i = . Y{t) (102) 



yii(i) = 30*2 + 9*^ + 0(t4) 
1 8 

Y2iit) = 360t - + 0{t^) 

Y31 (t) = -30*2 + + o{t^) 

010 

Y4i{t) = -360* + 180*2 - — *3 + 0(*4) 



1 R 

ri2(*)= 360*-yi=^+O(*^) 

F22 W = 4320 - 1080* + ^*2 - y*3 + 0(*4) 

y32(*) = -360* + 180*2 - ^*3 + 0(*4) 

8520 , 2148 o . 
1^42 (*) = -4320 + 3240* —t^ + -^t^ + 0{t^) 
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Y23{t) = -360i + 180^2 - H5i3 + 0(t4) 

>33(i) = 30*2 - 21t3 + 0(i4) 
1 242 

318 



8520 



2148 



Yu{t) = -360* + 180*2 - — *3 + 0(*^) 
y24(*) = -4320 + 3240*- 

l34(*) 

y44(i) =4320- 5400* 



*3 + 0(*4) 



1242 

360* - 360*2 _|_ ^^^3 _^ Q/^4\ 



7 

23640 



9852 



Note 



lim Y(t) 



0" 
4320 -4320 


-4320 4320 



(103) 



Now the original integral can be estimated by estimating the diagonalization 
of y(*). Though the proof docs not depend on knowing the origin of the t/(*) 
and -D(*) used to approximate diagonalizing !'(*), their construction is given 
in the appendix of the arxiv version of this paper [1] which also includes 
associated maple code and output. 



J(*) := 



776 det A{t) det A{t) 
1 



/ 



776 det A{t) det ^(*) J^a 



TT^det A(*)detyl(*) 



\x\\ 



lfe|2 



|x2|2|e 



-1/12 



|a;2|2| e" 



-t-^{Y(t)v,v) 

{Y(t){t-U),(t-iv))^^ 



Making the substitution w = t is actually saying 



Wl 
W2 
W3 
W4 

'dw\ 

dW2 

dw3 

_dW4_ 



*"2 



*~5 



'dhi 
dxi 
dh2 

dX2j 



dw = -dwidwi . . . -dw^dwi = * ^^-dhidhi . 
2 2 2 



-dX4dX4 = t 

2 
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making 



J{t) = 

1 



7r6 det A{t) det A{t) J^^ 
1 



30 



7r6 det A{t) det Jc* 



7c* 



Now 



Y{t) = U{tyD{t)U{t) 



4x4 



(diag) 

= U{trD{t)U{t)+[0{t')]^ 



J4x4 



where 



D{t): 



Ai 0' 
A2 

A3 
.0 A4. 

Ai 0' 
A2 
A3 

.0 A4. 



'0{t^^) " 

C(ii2) 


0(ii2) 



[0(i")]4x4 



(diag) 
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where 



T o..^ ^.o^ 25020 2 10050 , 66380 , 261767 , 
A. :=8640 - 6480 * + ^ - ^ ^3 + __ ,4 _ __ ,5 

48960935 g 29628553 7 208429618963 § 
2173248 * 8149680 ^ ^ 427173626880^ 
560822276587 g 46335059891 10 

~ 8543472537600 * 6133775155200 * 

518190034231 

1794129232896000 ^ 

r . s o 4 111 >s 161 fi 20561 7 561019 . 

A2 :=6t^-3rH f V ^ r 

80 960 1209600 21772800 

3916753 9 _ 827998967 

15676416000 * ~ 282175488000 * 

5185091420987 ^ 

15643809054720000 * 
~ 141.1.17 6223 c 256685 q 

A3 :=- r e + — f ^ — r H e h r 

3 12 72 32 207360 8957952 
588107563 ^ 6399891227 „ 

22574039040 325066162176 

A.:=^t^-lt^-^.^-i^.«^^.^ 
8 16 768 3072 40960 

37303 10 _ 2563021 
~ 1474560 ^ ^ 123863040 * 



and 



where 



[/(O = U{t) 



oitA 



o(A 



om om 
om om 
om om 
om o(A 
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V2 V2. 
Uii = 1 1 

24 48 

1 3 65 2 
Uai = 1 1 

2 16 2304 

V2 „ 43^2 2 

13 193 2 
Uii =- — —t — 



16 



768 



29 



\/2 \/2 5\/2 , 

Ui2= - —t + — 

2 8 288 

13 5 o 

= 1-\ r 

^^ 2 16 256 

„ 6IV2 2 
11 53 2 

Ua2 = 1 1 

2 16 768 

11 167 2 
^^^=-2 + T6*+ 2304* 
V2 „ 7^2 2 

11 215 2 
^^^=2 + 16*+ 768* 

V2 \/2^ 5\^^2 

Wl4 =— —t TTTrrrt 

2 8 288 

11 41 2 

W24 = — I 1 1 

2 16 768 

V2 109V2 2 

"^^~*+^* 

1 3 29 2 

=2 + 16* - 768* 

and taken so that U{t) is invertible. 

The following lemmas say that U{t) is approximately orthogonal. 

Lemma 2 U{t)U{tr =1+ [0{t')]^^^ 

Proof 

U{t)U{t)* 

4x4 

= [mt)+[o{t')],^,)(uitr+[o{t^)]:^^^ 



30 



where 

^ ^ 1152 13824 
' 9216 36864 

^ ^ 27648 110592 

4,. / X , 13 3 12317 , 

M22(t) = IH H 

^ ' 2304 1327104 

^ 1024 165888 

^ ^ 1152 18432 

^ ^ 1152 13824 

' 1024 165888 
4M33m = l + — + 

' 192 4608 

' 9216 110592 

^ ' 9216 36864 

^ ^ 1152 18432 

4M34(t) = ^ 72*3 _ ^^4 

9216 110592 
4,. / N . 95 o 21781 4 
^ ' 768 147456 

so 

Lemma 3 f/(t)-i = f/(t)* + [©(t^)]^^^ 
Proo/ 

^ f/(tr = [/(r^(7+[o(i3)],^j 
^ f/(r^ = c/(tr(/+[o(f^)],^J"' 



□ 
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since [0[t-^)\^^^ > V^iixA 



so 



u{t)-' = uitr + u{tr[o{t')]^^^ 



□ 



Now 
J(i) = 



f 



[ ||t«i|'-Hn-|k3|'-Kne-<^W*^«^««''"'>d«; 
7r6 det ^(t) det ' i IM i | 



TT^ det^(t) detyl(t) 

^30 



Make the substitution 



z 


:= U{t)w 






Wi 


3 


dw 


= det{U{t)) 



dz + 0{t^) dz 

In the following, only the properties of the t^, and terms of the 

used so, for the sake of readability, "w^i" will always be written in place of 



J{t) = 

^30 



77^ det A{t) det A{t) 



(105) 
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+30 



TT^ det A{t) det A{t) Jc* 



33\ 



TT^ det A{t) det A{t) Jc^. 



33\ 



TT^ det A{t) det A{t) Jc* 



1 V- 


2 




2 1 


1 z2j ujizj 








ll „ 


2 




2 


1 ^3 




N , 7/ z^' 






2 












2 1 




2 




2 






"34 






2 












2 1 




2 




2 






Ej «34 2j 





-(n(t)z,z> 



e ^^=1 '^^^^^^ d2; 



Make the substitution 



'^3 '■= V>^Zj 



and dwj = 



so 



1 _ i _ i dwi dwi i dwi dw^ 
dz = -dzi A dz\ A ... A -dz^ A dzi = - A A ... A - A —i= 

2 2 2 yAi V Ai 2 y\i V A4 

1 i z 
= -dwi A dw)i A ... A -cJw4 A dwi 

A1A2A3A4 2 2 

dw 



A1A2A3A4 



so, since Lemma 2 and (100) imply 



[detyl(t)]AiA2A3A4 = t^" + ©(t^^) 



J{t) = 

^30 + 0(^33) 
7r^detA(t)detA(t)nA 



^30 + 0(^33) 



7r^detA[t™+0(t^^)] Jc* 
^10 + 0(^13) 



■ ||s.-"i3^|-|E,-«.-4^|'|e" 

^ . I 9 I _ . I 9 I 



I ^."^-7^7 1 



7r^detA(l+0(t^)) 



i + o(t3; 



• ||E,«.3^|'-|Ei«.4^|'je" 



TT^ det A{t) Jc4 
l + 0(<3) 



* ||Ej«ji^| -|Ei«.-2^| 

.51 



■K^ det A{t) Jc4 



I 

Jc 



t2 , 



\A7 



I ^3 ,A 



Now 



,5 

f2tm 

5 

t2?il2 



VA7 



30 7 7V30 9 ^/ 11 



8640 
V30 



720 

30 



i2 - 



69120 
/30 7 

t2 



VAT 

r2 ■U]^4 

AT 



8640 
^30, 



t3 



1152 

30 



720 



i2 



23040 

V30 
5760 



+ o(ii) 



34 



and 



V6 
12 



12 



^2^31 



\A3 



v6 1 v6 3 

— i2 + —t^ 

2 16 

•\/6,a 67-\/6,4 



576 



+ o(ti) 



— ^ = -^t^ - -Tzrt^ + Ou2 
VA^ 2 16 \ 

t^U34 3 121-\/6 5 



V6 

6 ' ^ 576 



so 



J(i) = 

i + o(t 



^ V6 _ 7V6 473V6 2 (^3) 
y/Xi 3 72 3456 ^ ^ 

^^M42 _ 29^6 2 Q/,3\ 

VAI 3 72 ' 3456 ^ ^ 

^5^43 5-\/6 799-\/6 9 ^/ ox 

= — + -^i H + O(t^) 

VAI 3 72 3456 ^ ^ 

f5u44 \/6 ll\/6 91-\/6 9 ^/ 



• hitiwhuiw) - 6u{w)6u{w)\ ■ e-<"''«'> dw 



35 

where 

ait{w) = 0{t^)wi 

+ ( + 7^*^ +^(*^) ) ^3 

fVQ 7V6 473^6 2 ^/ 
fVQ Ve 29V6 2 ^/ -hnA 

(•\/6 1 -\/6 3 ^, 5 , \ 

/V6 5^6 799^6 o ^/ 
(5it(w;) = 0{t^)wi 

fVe ll^/6 91\/6 2 ^/ s\A 
V 3 72 3456 ^ ^ ^ 
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So 



jk 



WjWk 



Ui + /34\/i + 74* + ^4*2 + ^^£jfc {t)wjWk 

jk 



where 



a2(i) = — ■u;4— w;4 —Wi—Wi = 



_ , , VQ \/6 _ \/6 \/6 _ , _ _ N „ / - X 
P2(i) = ^'^3-^W4 + -^^4^ty3 = + mm) = 2Re(u;3W;4) 



^3(i) = ^^3-^W4 + 



W3——W4 = — (W3W4 + W3W4) 



= -2Re(w3w;4) = -/32(t) 



, \/6 , 7x/6,_ 



7^6 VQ 

W4- 

V6 



72 ^^l""'^ 



72 



W4 



Vg Vq _ 



= - (— 2W2«'4 — 2w2Wi — 2WiWi + 9W3W3) 

= ^(-4Re(w;2W)4) + ^\m\'^ - 2\wi\'^) 
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\/6 \/6 _ , VG, I \/6 _ , a/6 V6, _ 

= - ^^2-^W4 + ( ^)W3( ^jU'S + — W4(- — )u;2 

\/6 5\/6_ 5\/6 y6_ V6 \/6 _ 

\/6 a/6_ V6 11V6_ ll\/6 V6_ 

;;-W'4-7Tr'"'2 — W4 „^ W4 ^'4— — W4 

312 3 72 72 3 

= - (—2^2 ^^4 — 2«)2W4 — 2t«4 tl;4 + 9^3 tDs) 

=i(-4Re(u;2i(>4) + 9|w;3|2 - 2|M;4n 
= 72(0 



1 _ _ _ _ 3 

<52(i) = -Y^(3u'2l«3 + 3W2W3 + '^W-iWi^ + 4ui3Ui4)t2 

= -^(6Re(t«2'it>3) + 8Re(w3«)4)) 



1 _ _ , 

= ■^{6Re{w2W3) - 16Re(w;3it)4)) 



notice 0:2, Q!3, ,02, /Ss, 72, 73, ^2, (^3 G K. This implies that since 
^a2 + /32\/i + 72* + ^2*^ + ^^ejfe(i)wj«)fe^ e 11 



and 



Q!3 + ^3\/i + 73* + ^3*2 + y^\jk{t)WjWh 



jk 



by (105), (^jk^jk{t)wjWk^ and (^jk^£jk{t)wjWk^ must be real, as well. 



38 



so 



J{t) = 

l + 0{t^) 
7r6 det A{t) 



jk 

a4{w) < I34iw) 



TT^ det A Jo^i 



jk jk 



where 



a4 = -(^2)^ = -{2Re{w3W4)f = -4Re(w3W4)^ 
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/34 = /32<53 + 72^ - S2I32 

= — Ke{w3W4) ^6 Re(u'2tZ)3) — 16 Ke{w3W4) 
+ ^ - 4Re{w2W4) + 9\W3\'^ - 2\w4\^^ 
+ - Re{wsW4) ^6 Re{'W2W3) + 8 Re{wsW4) 

= 2Re{'W3W4)Re{w2W3) - ^(Rc(w3W4))^ + ^{Re{w2W4))'^ 
4 9 o 2 

— 2 Re{W2'W4)W3iB3 + - Re{W2'W4)'W4iU4 + -Ws iV^ — Ws'W3'W4W4 

1 o _o 
+ qW4 W4 

= 2Re{w2Wz) Re{w3W4) — ^ Re(w3W4)^ 

+ ^ Re(w2W^4) ^4Re(w2W'4) — ISl^sl^ + 4|'u;4| 

+ l\ws\'-\Ws\-'\w4f + l\W4\' 

= 2 Re{w2W3) Re(u'3W4) — - Re(w3w)4)^ 

+ ^ Re(w2W'4) ^4 Re(w2'it'4) — 18|w3|^ + 4|w4|^ 

+ -^(9|«;3r-2|«;4|') 



3 

4 



Again, notice that a4(w), /34(w), i73(w,t), and S4{w,t) are all real. 
So 

J{t)= It? 1*11 / |a4t + /34i' + i:3K0 + ^4Kt)| e-<'"''">du; 

The following lemma will help take the error terms out of the absolute 

vahic. 

Lemma 4 For any n > 1, any continuous complex valued f{w,t) and g{w), 
and s{t) = 0(i") 

\f{w,t)\ = \f{w,t)-eit)g{w)\+0{n\g{w)\ 
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Proof 



\f{w,t)\-\f{w,t)-e{t)g{w)\ 



< 



e{t)g{w) 



m\ \9{w)\ 



so 



H£(i)| • \9{w)\ < \f{w,t)\ - \f{w,t) - e{t)g{w)\ < \e(M ■ IflHI 
0(i") 0(i") 



so 



so 



\f{w,t)\ - \f{w,t) - e{t)g{w)\ = 0{n \g{w)\ 



\f{w,t)\ = \f{w, t) - e{t)g{w)\ + 0(r) \g{w)\ 



□ 



Applying Lemma 4 repeatedly to \a4,t + (iit^ + Si{w,t) + Si{w,t)\ and 
\jk{t)wjWk says that 



J{t) = 

l + Oif) 
7r6 det A{t) 



and applying Lemma 4 repeatedly to \a4t + fiit^ + Ei{w,t)\ and 
^s:jMm{t)wjiBkWeWm says 



jkim 
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Finally, applying Lemma 4 to \a4t + /34<^| and /34t^ says 

7rMetA(t) Jci \ ^ ^ ' 

JK 

jkim ' 



so 



l + 0{t^) ( 



7r6 det A{t) 



i + o{t^) 

7r2det^(t) 



( -t(-27r*) + 0(t2) [finite] [^^'^e 
+ ^ O(i^) [finite] 



jkim 



{2t + 0{e)+0(ti)+0{t^)) 



^2dctA(t) + 
2t + 0(t2) + 0(t3) + + 0(^5) 

7r2det>l(f) 

2i+0(t2) 
7r2detA(t) 
1 2t + 0(t^) 
^ 3t + 0(i2) 

= 3^ + ^(*) 



37r2 
i.e. 

J(r) > 



37r2 



8 Appendix 

or Where did U(t) and D(t) come from? 
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D{t) and U {t) naturally arise when calculating J{t) by diagonalizing Y{t) . 
1 



TT^ det A{t) det A{t) J^a. 



I / 

776 det A{t) det A{t) J^a 

776 det A{t) det yl(t) ic* 



/ \\h,\^-\x,\'\.\\h2\'-\x,\^\ 

.g-{y(t)(t-i^),(t-5^)>^^ 



776 det A{t) det Jc* 
776 det det A{t) Jc^ 



\\Wl\ - \W2\ I • \ \Ws\ - \Wi\ I 

• e- dw; 

. g-((7(*)*_D(t)C/(t)™,,«) 



\WI? - \W2?\ ■ \\Wi? - \Wi?\ 



-{D(t)U(t)w,U(t)w) 



where 



D{t) 



Xi{t) ■ 

X2{t) 
Xsit) 
X4{t) 



with the Xi{t) being the eigenvalues of Y{t) and 



U{t) 



V2{t) 
V3{t) 
Viit) 



with the Vi{t) being the associated normalized eigenvectors, making U{t) real 
orthogonal. 

So expanding the Aj and Vi in t gives an expansion for J{t). Maple out- 
puts a 100 megabyte file for each eigenvalue when asked eigenvalues(Y(t)) 
directly and crashes when asked to find an expansion for any individual eigen- 
value. However, because the matrix is 4 x 4, the eigenvalues and eigenvectors 
can be calculated algebraically by applying the quartic formula [7] to Y{tys 
characteristic polynomial. Then these algebraic expressions can be expanded 
by maple. 



Y{t) =: 



7iW /2W /aW fA{t)- 

Mt) hit) frit) fsit) 

hit) hoit) hiit) h2it) 

./isW hiit) h^it) /leW. 



43 



det(y(t) - xl) = 











-fifefie^ 




^ — /sA/ii/ie"^ 










+/13/6/4 




— /5/2/12/15 










"/i/e/ii 




— /5/10/3/16 










— /5/14/4 




+/5/10/4/15 










+/1/10/7 




+/5/14/3/12 










— /13/3/12 




— /5/14/4/11 






/ -/9/3\ 




+/13/4/11 




— /9/2/7/16 






+ /6/16 




+/5/2/11 




+/9/2/8/15 






— /13/4 




—79/2/7 




+/9/6/3/16 






+/6/1I 




+/9/3/16 




— /9/6/4/15 


/ -/6\ 




-hf2 




— hfiflh 




— /9/14/3/8 


-A 












+/9/14/4/7 


-hi 


+/11/16 


+ .f5.f2.fl6 


+/i 3/2/7/1 2 






— /12/15 




— /1/11/16 




— /13/2/8/11 






— /10/7 




+ /1/12/15 




— /13/6/3/12 




+hh 




— /13/2/8 

+79/6/3 




+/13/6/4/11 

+/13/10/3/8 






\ +/1/16/ 




+ /1/14/8 




—/i 3/1 0/4/7 










— /e/ii/ie 




— /i/e/ii/ie 








+/6/12/15 
+/10/7/16 

— /10/8/15 
-/14/7/12 




+/1/6/12/15 
+/1/10/7/16 

— /1/10/8/15 
— /1/14/7/12 










\ +/14/8/11 / 




V+/1/14/8/11/ 



Because F is hcrmitian, it's eigenvalues will all be real. In fact, because K 
is positive definite, they will all be positive. In the calculation below, things 
can be complex (e.g. i?) but all of the imaginary parts will go away by the 
end. 

F3(i) = -8640 + 0(t) 
F2(i) = 51840^3 + Ci(t'*) 

= -17280i^ + Ci(t*) 
Fo(i) =6480t^^ + C>(t") 
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Solving the general quartic a;^ + Bx^ + Cx^ + Dx + E = requires some 
simplifying definitions and a few choices. 

a := — +C = -27993600 + 0{t) (< for small tast\0) 

8 

/3:=- —+D = -80621568000 + 0{t) (< for small i as t \ 0) 

8 2 

CB^ BD 

7 := -— + — — + E = -65303470080000+ 0(t) 

Q,2 

P := -— - 7 = -223948800^*^ + 0{t^) 

Q:=--^ + ^- ^ = - 1289945088000 + 0(t^°) 
108 3 8 



^ = -2 ±VT+27 

= 644972544000*'' + 0{t^°) ± ^-8666449635704832000000*20 + 0(*2i) 
(Choose either root.) 



iO(tio) 



[/ — (Choose any of the three roots.) 

5 rr P 



W := ^Ja + 2y (Both roots come up in zt^) 

With those definitions, the roots should be 



±sW±t J- (3a + 2y±, f) 



B 

-4+ 2 

where the ig's are dependent and the ±t is independent. 

Maple can expand a, /3, and B easily. We need to carefully intervene to 
get it to expand W, y, and ^ which require R and U. 

~ ~i ~ 27 ~ 

on 563370492281772061551028318944792605312337 An ^ / 41 \ 

8666449635704832000000 t -\ H 1 +0{t ) 

7662929083743683750 ^ ^ 
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g2 p3 



A/86664496.'i'')704832()()0000t2o h ^ 0(t4i) 

27 



=\/8666449635704832000000 i 1° • ^1 + • • • + 

^ V 



=\/8666449635704832000000 f ^° ■ (l + ^x + ---+ 0(a;^^) 



=\/8666449635704832000000f^° + • • • + 0{t^'^) 



= 644972544000 1'^ + • • • + 0{t^°) 



El 

27 



2 



V8666449635704832000000ti° + • • • + 0(^31) 
644972544000 19 + • • • + O(<30) 
1 1 



644972544000*9 H + 0{t'^'^) 



( -1/8666449635704832000000 + • • • + 0{l?^) ) 



1 

l+x 



-(l-a; + --- + 0(x2i)) 



644972544000 

• (\/8666449635704832000000 + • • • + 0(t'^^) ) 

: \ il + ■■■ + Oit^^\ \ 

644972544000 t9 ^ ^ 

■ (\/8666449635704832000000 + • • • + 0{t^^) ) 



= 644972544000^9 ( V8666449635704832000000t^° + • • • + 0{t'') ) 
Choose R in the I quadrant, namely 



2 V 4 27 
=644972544000 + • • • + 0{t^°) 

+ i(\/8666449635704832000000i^° + • • • + 0{t^^) ) 
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where 



/_p3 

tr = \l and Br = arctan ^ ^ 



Choose U in the I quadrant, namely 
where 



Q2 _ p3 

4 27 



7[/ = - arctan -^^ and w ~ 



'2 



^ 6 3U 6 Srije'^c; 



Pe-iOu 

-tQ! + me""^ 

6 3r(7 

-«"+-^ 

6 ru 

5 rfj{cos6u + ism6u) — ^{cos9u — isinOu) 
--a-\ ^ 

6 ru 



6 ru 



H + f ) 



ru 



727 



\ 2 



COS 0u + i 



sin 

2 



ru 



sin 0[7 



.5 

6 \ ru 



^ cos 6'i7 + i 



^ ^ sin^t/ 



ru 



5 2P 

6 3r(7 



cos 6u + i \ I sin ( 



5 

— rv 



2P 
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so 



arctan ^ — ^q"^ j = arctan H h 0(t^^)^ 



larctan(VZ^).v|, + ... + On 



2 



COS - arctan I ^ 1 = cos y-^t H 1" 0[f'^)j 



2 

l + ...+0(x22) = l + ...+0(t22) 



= - ^ (-27993600 + • • • + O (f^-^) ) 
6 



+ —= (ssAoVst^ + ...0(t2^)) (l + . . . + o{t^'^)) 



5.27993600 ^/ 2>s\ 2 • 8640\/3 , ^/ 2>;\ 

= ^ + ... + o{t'') + — -^t' + ... + o{t'') 

5-27993600 

= + ... + 0(f25) 
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= ^(-27993600 + • • • + 0{i^^)) + 2 + . . . + 0(f^) 

=\ - ■ 27993600+ • • • + 0(t25) = W- • 27993600^/1 + • • • + 0(i^^) 
V 3 V 3 ^ ' 

X 

• 27993600 (! + ■■■ + ©(.t^^^)) = • 27993600 (! + ••• + 0(i^^)) 



'I • 27993600 + • • • + 0(i^^) = 4320 + • • • + 0{f^) 



2/3 -161243136000 H \-0(t'^^) 

W ^ 4320 + • • • + OCt^s) 

= (-161243136000 H h 0{f^)) 



4320 ! + •••+ 0(i25) 



1 

l+x 



-(-161243136000+ • • • + 0{f^)){l + ■■■ + 0{x'^^)) 



4320^ 

161243136000 o.x ^/ o^x 
+ • • • + 0(t^^) = -37324800+ • • • + 0(t^^) 

4320 ^ ^ V y ^ ^ V ; 
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So 



B W + ^-{3a + 2y+§) _ -B + 2^ -3a - 2y - § 



17280+ ••• + 0(t25) + 17280(1 + ••• + 0(a;25)) 



4 



8640 + ••• + 0(i^^) 



^^^=^-■+ = -4+ 2 4 

/ -(-8640+ •••+0(^25)) + N 

= "7 / — 3(-27993600H i-o(t^-')) 

^ +2W ^ ^ 

y y - 2( 5-27993600 ^...^o(t25-)) „ ( - 161243136000 ^ ... ^^(^^25 ) ) J 

1 /(8640 + -- • + 0(^25)) +2(4320 + - •• + Ci(t25)) \ 
4 1^ + 2^74649600 + • • • + 0{t^^) J 
_ 17280 + • • • + 0(i25) + 2 • 8640v'l + • • • + 0{t^^) 



-B + 2W -2^-3a-2y-^ 
■- A+,- ^ 

_ 17280 + ••• + 0(t25) - 17280(1 + ••• + 0(i25)) 
~ 4 

= Qt^ + ■■■ + 0{f''') 



A3 := A_ + = 



-B-2W + 2\l-3a-2y+^ 



-B-2W + 2y^it8 + --- + 0(f25) 
4 

-B - 2T4^ + ftVl + • • • + <5(*^^) 



-B-2W^+-i4(l + .--+0(xi^)) 



-B-2W+lt^ + -- 



o{i 



21\ 
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A4 := A_,_ = 



-B-2W 



-3a-2y+f 



_B-2W-|t''(l+-+0(fi^)) 



We can check that these expansions for the eigenvalues are correct by 
checking that evaluating the elementary symmetric polynomials in four vari- 
ables on them gives the coefficients of the characteristic polynomial, i.e. 



characteristic polynomial = 

- ei(Ai, A2, A3, A4) + e2(Ai, A2, A3, X4,)x'^ 



e3(Ai, A2, A3, A4)a; + e4(Ai, A2, A3, A4) 

^ V ' 

det 



Maple confirms that 



ei(Ai,A2,A.3,A4) + B = 0(P) 
e2(Ai,A2,A3,A4)-C = 0(i2i) 
e3(Ai,A2,A3,A4)+i?-0(i24) 
e4(Ai,A2,A3,A4)-i; = 0(t2«) 



To find the associated eigenvectors, we need to find the kernel of A := 
Y — Xjl. i.e. we want to find Vj such that Avj = 0. If the last row of A isn't 
a row of zeroes, which is the case for the four A's that we examine here, then 
a sequence of elementary row operations represented by multiplication by an 
invertible J can take ^ to a matrix with a row of O's on the bottom. 



JA = 



B 
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then 



JA 



-B- 



-BB- 



B 



a 
b 
c 







a 




a 






"0" 


1 


b 


+ 


b 











c 




c 


























-B- 





a 






"0" 




"0" 




b 




= J-i 












c 












1 

















a 
b 
c 



so 



-B- 



is an unnormalized eigenvector. We need —B ^ 



for 



each A. To expand B ^ we just expand adj B and det B and divide to get 



^.1 vi: 



A = 
-8640 + • 
360t + -' 
-30t'^ + - 
-360i + - 
We get vi = 



• + 0(i^^) 360t+-+O(t^^) -30t=+-+O(t^=) -360t+-+O(t^=)' 

■+0(t^^) -4320+---+O(t^^) -360t+---+O(t^^) -4320+---+O(t^=) 

■ + 0{t'^^) -360t+---+O(t^^) -8640+---+O(t^=) 360 t+---+0(t^^) 

■ • + O (i^^) -4320+---+O(t^^) 360t+---+O(t^^) -4320+---+O(t^^) 



■_\/2 , 
24 



72 .2 I 

48 ^ 



V2 _ V2-I. I 
38 8 
# + 0^2 + . 
\/2 _ \/2 



+ 0(t25) 



Maple confirms that 



Yvi — XiVi 



'0 






o 













o 


) 
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8.2 V2: 



A = 



30<^H hO(t^-'^) 360 «H \-0{t^^) -30<^H hO(t^^) -360 tH |-0(t^'')' 

360 tH |-0(t^'') 4320+---+O(i^''') -360 t+---+0(t^-°^) -4320H \-0{l^'') 

-SOt^A hO(t^^) -360tH hO(t^^) 30t^ + ---+O(t^-'') 360 t+---+0(t^^) 

-360 tH hO(t^^) -4320H |-0(t^'^) 360 t+---+0(t^^) 4320H |-0(t^^) _ 



We get V2 



2 16 ^ 
1 3 



■+0(tl8)' 

+ o(ti«) 



t + - 



16. 



again, maple confirms that 

Yv2 - X2V2 



8.3 V3: 

"SOt^H |-0(t^^) 360t+---+O(t^^) -SOt^H |-0(t^^) -360tH |-0(i^^)' 

360t+---+O(t^^) 4320H |-0(t^^) -360tH |-0(t^^) -4320H \-0{t^^) 

-SOt^H l-O(t^i) -360tH l-O(t^i) SOt^H \-0(t^^) 360t-\ \-0[t^^) 

_-360tH |-0(t^^) -4320H \-0{t") 360tH |-0(t^^) 4320H \-0{t") _ 

In this case, tlic unnormalizccl eigenvector would have negative powers of t, 
so in the maple calculations, we find the normalized version of tv which is 
the same. i.e. we use the fact that 

tv V 

We get t;3 = 

^^ + o^ + ... + o(^9) 
^t + ^^t^ + ... + o{t^) 
_^^ + o^ + ••• + o(^9) 

again maple confirms that 



Yvs - A3V3 = 



"0 






















_ 
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8.4 vr. 



A = 



30<^H hO(t^^) 360 <H |-0(t^^) -30<^H hO(t^^) -360 tH l-O(t^i) 

360tH |-0(t^^) 4320+---+O(i^^) -360 t+---+0(t^^) -4320H 1-0(<^^) 

-30t^H hO(t^^) -360tH [-0(t^^) 30 t^ + ---+0(t^^ ) 360 t+---+0(t^^ ) 

-360tH hO(t^^) -4320H |-0(t^^) 360 t+---+0(t^^ ) 4320H |-0(t^^) _ 



We get V4 



1 

! 

2 
1 

2 
1 

L2 



1 

16 



2 ' 16" ' 



+ 






+ 






+ 


r< 




+ 







maple confirms that 



YV4 — X4V4 

\ 1 collected together, we have 



Ai 


= 8640 + • 


• • + o{t^') 


A2 


= 6t^ + • • 


• + 0(i25) 


As 






A4 







and [/ = [[/ij] where 



f/21 = -i-^.+ - + Or) 

C/3i = ^ + 0i + .-. + O(f9) 
U4i = l:-^t + --- + 0{t'') 



U22 = 
U32 = 



V2 6IV2 2 ^/ 



Ul3 
U23 



2 16 
V2 



?742 = 



1 



1 



24 
1 



t + Ot'^ + ---+0{t'^^) Ui4 

1 



16 



t + --- + o{t 



18\ 



U24 



V2 



U33 = -^ + 0t + --- + O{f) U34 



2 16 

~2 8 

1 1 

2+16* 

V2 109V2 

—t H — 

6 576 



+ o{t 



18\ 



16 
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The calculation only requires the Aj be expanded to t and U be expanded 
to t^. 



Acknowledgements I would \ike to thank my advisor Bernard Shiftman for his 
advice, patience, encouragement, and many ideas about solving this problem. 
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